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ABSTRACT 


Integrated and spectral error energetics of the GLAS General 
circulation model are compared with observations for periods in January 
1975, 1976, and 1977. For two cases the model shows significant skill 
in predicting integrated energetics quantities out to two weeks, and for 
all three cases, the integrated monthly mean energetics show qualitative 
improvements over previous versions of the model in eddy kinetic energy 
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and barotropic conversions. Fundamental difficulties remain with 
leakage of energy to the stratospheric level, particularly above strong 
initial jet streams associated in part with regions of steep terrain. 

The spectral error growth study represents the first comparison of 
general circulation model spectral energetics predictions with the 
corresponding observational spectra on a day by day basis. The major 
conclusion is that eddy kinetic energy can be correct while significant 
errors occur in the kinetic energy of wavenumber 3. Both the model and 
observations show evidence of single wavenumber dominance in eddy 
kinetic energy and the correlation of spectral kinetic and potential 


energy 


preceding pa ge 
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1 . Introduction 

Analysis of the generation, transport, and dissipation of energy 
represents one of th* primary methods of studying geophysical fluid . 
flows. Lorens ( 1 955 ) proposed subdividing atmospheric energy into 
kinetic and available potential forms and in turn subdividing these into 
zonal and eddy categories. Subsequent observational (Oort, 1964) and 
general circulation sodel (Kanabe and Terpstra, 1974; Kasahara and 
Washington, 1971; and Somerville et al., 1974) studies have followed 
Lorenz's approach and will be referred to as integrated energetics 
analyses. Further subdivision in the spectral domain was suggested by 
Saltzman (1970) and followed by similar observational and model studies 
(Tenenbaura,1976; Wellck et al., 1971; Baker et al.,1978). One purpose 
of this paper is to examine both the integrated and spectral energetics 
of the current Goddard Laboratory for Atmospheric Sciences (GLAS) model. 

The most direct test of a model over a period of several days is 
the quality of the prediction of standard meteorological variables. 
Comparisons are made either in terms of synoptic discussions (generation 
and deepening of lows, etc.), root-mean-square errors, or the related 
idea of skill scores (Tewles and Wobus, 1954; Atlas, 1979). One knows 
theoretically that these approaches should give an essentially random 
result at a predictability limit of the order of 14 d and in practice do 
so in about 5 d (Lorenz, 1967). 
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To retract further information in spite of these limits, one 
employs various space and time averages applied to both meteorological 
variables and energetics quantities. These averaged quantities may show 
skill for periods approaching or exceeding the predictability limit; 
indeed current general circulation models are characterised by quite 
reasonable model "climates. " One of the commonly held intuitions is that 
low wavenumber spectral quantities qualify as "averaged" quantities. A 
second purpose of this paper is to examine that assumption and to 
display circumstances where it is more or less valid. 

In this paper we will follow Saltsman in employing one-dimensional 
wavenumber spectra consisting of Fourier transforms around latitude 
circles. An alternate analysis in terms of two-dimensional spatial 
spectra was suggested by Baer (1972) and is related to many spectral 
models. A subsequent paper will present two-dimensional results 
(Tenenbaum, to be published). 

The wavenumber domain spectra presented here form one member of a 
trio of approaches concerning atmospheric wave y '.enomena. Frequency 
domain studies deal with spectral analyses of variable* in that domain 
while space-time analyses take both approaches simultaneously. These 
approaches are complementary, with the wavenumber approach allowing us 
to study the time dependence of energetics and spectral quantities. In 
particular, we can study the growth of errors toward the predictability 
limits for a variety of synoptic situations. It is the study of these 
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synoptic dependences which form ths third purpose of this paper. 

In this paper we do not separate transient and stationary effects. 
As we shall subsequently see the intense stationary jet stream anchored 
east of Japan tends to dominate wavenumber 3 results. While time means 
could be subtracted from each of our independent variables, we ure 
hesistant to do so we have a better sense of the interannual 

variability. A final ; nd distinctive feature of this study is the use 
of model spectra in comparison with the spectra derived from obser- 
vations for corresponding times. Studies of model spectra and of 
observational spectra have been separately performed; no spectral 
energetics study for the corresponding times has yet been published. 


2. Model end methods 


The GLAS general circulation model represents the third in a aeries 
of mo***'*. The original Mints-Arakava model (Arakava, 1972) was a 
3-leve' todel incorporating Arakava’ a energy and enatrophy conserving 
difference scheme. The second in the series, the Goddard Institute for 
Space Studies (GISS) model vaa a combination of Arakava* s approach vith 
modified convection and radiation echemes (Somerville et al., 1974). 
Subsequent papers have dealt vith its spectral properties (Tenenbaum, 
1976), synoptic performance (Druyan et al., 1974), seasonal behavior 
(Stone et al., 1977), and behavior over aemi-arid regions (Charney et 
al., 1977). 

The CLAS model differs from the GISS model primarily in the 

application of a Shapiro filter to the dependent variables and an 

improved radiation routine. Specific changes include (Halem et al., 

1978): (1) improved l*ng vave radiation formulation, (2) increased 

radiative process time step, (5) prognostic rather than prescribed soil 

moisture along vith other hydrological changes, (4) smoothly varying 

sea-surface temperature changes, (5) changed surface albedo, (6) 

coarsening of tonal resolution in 5 bands tovards the poles, and 

% 

ine. eased advective time step (10 min) vith decreased smoothing near the 
poles. 


The integrated energetics quantities are calculated according to 


the formulas given in Oort (1964) for the mean tonally averaged kinetic 
energy, K M , and the tonally averaged eddy kinetic energy, Kg (l.e. the 
tonal average of the deviations from the tonal mean). Using Oort* a 
nomenclature, both quantities are calculated in the space domain. 

The spectral analysis methods are taken from Saltsman (1970) and 
are described in further detail in Tenenbaum (1976). One-dimensional 
Fourier spectra are calculated around each latitude circle at each of 
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the nine levels every 12 h. Tropospheric results represent the average 
of the lowest 8 levels (centers from 175 mb to 945 mb) over the entire 
Northern Hemisphere. Because of problems with National Meteorological 
Center (NMC) observational results at 1200 CMT prior to April 1975 (Spar 
et al., 1976), all graphs are from 0000 GMT data. No significant 
differences appeared when comparisons also Included valid 1200 GMT data. 
The three periods dincusaed in thin paper represent model predictions 
initialised at 0000 GMT 1 January 1975, 1976, and 1977. The climate 
averages represent 1 month averages from 0000 GMT 1 January through 1200 
GMT on 31 January. The time dependent graphs run over the periods 
noted. 


3. Synoptic energetics summary 

Ve present in this section a synoptic summary of the periods 
covered, in both conventional and energetics terms. Tropospheric 
a results refer to the Northern Hemisphere troposphere, consisting of 
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levels 2 through 9 of the CLAS model and extending from 120 mb to the 
surface. Integrated energetics diagrams have been presented by Oort 
(1964) for observations and each of the major modeling groups for the 
general circulation models (Manabe et si., 1970; Vellck et al. f 19 7 1; 
Stone et al., 1977). The typical time dependent behavior (Baker et al., 
1978; see also Fig. 6 below) shows two primary features: Kp| and Kg 

tend to be anti-correlated and Kg typically varies between 9 and 9x10^ 
J/m^. 

4 

Spectral energotics distributions have been similarly documented 
(Tenenbaum, 1976) and show two characteristic features in the eddy 
kinetic energy spectrum. First, the bulk of the energy is in wave- 
numbers 1 through 9 and, second, there is a characteristic log-linear 
falloff for wavenumbers above 8. The relationship of the letter to 
theories of turbulence has been extensively discussed in the literature 
(Chamey, 1971). 

Since our subsequent concerns deal with the distribution and 
redistribution of eddy kinetic energy among the individual wave numbera, 
it is instructive to match the synoptic maps with the spectra. Fig. 1 
ahows the 200 mb map for 0000 GMT 1 January 1977. To a good approxi- 
mation the bulk of eddy kinetic energy is located at this level (Fig. 19 
in Tenenbaum, 1976 and subsequent figures below). Table 1 presents the 
corresponding energy apectra. 

Pxaminatlon of the map shows a very characteristic subtropical jet 
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•outh of Japan. Thia Jet apptara in conjunction with laaaar Jets 
eantarad ovar tha Western Atlantic and Saudi Arabia and ia conaiatant 
with tha quasi-stationary wintar average wava pattern found by 
Krishnamurti (i960. For January 1976 (not ahown) wa have a strong 
polar Jet eantarad ovar Norway and a such weaker subtropical jet in that 
aactor. These Jata provide tha dominant contributions to tropospheric 
Kg, and especially to tha spectral subcomponents of Kg, given by K n . 

Tha model's ability to predict tha growth and decay of these Jets will 
determine its ability to predict Kg and K Q . 

We can correlate the Jet stream behavior with the spectral analysis 

by comparing the 200 mb charts with the individual K given in Table 1. 

n 

For 1977 we see an initially dominant n ■ 2 pattern characterised by the 
intense Japan Jet matched with a North American Jet not quite 180° away. 
By 0000 GMT 7 January (Fig. 2) we see a flip to dominance by n • 1 and 
3 in Table 1 and a corresponding change in the Jet stream pattern. The 
Atlantic portion of the North American Jet has been cut off and a Saudi 
Arabian Jet has grown to above 150 KT (77 m/a). Two of the three peaks 
contributing to the n • 3 component are eaeily visible while the n • 1 
component arlsea from the leas obvious polar asymmetry. In more 
conventional tame the first half of 1977 wae characterised as follows 
(eee Wagner, 1977): A strong blocking ridge was located Just offweswm 

North America. This ridge weakened during the second week simultaneous 
with the growth a strong trough in t ' is Alantlc. This last effect 
occurred in conjunction with an Intense stratospheric warming whose 
effecte include a significant sonal flow over Hudson Bay (Fig. 3). 


9 


4. Integrated energetics 


Probably tha aoat coaaonly aeen fallura of general circulation 
aodala haa baan tha daficit of K fi . Moat aarly aodala had thia problem 
and whila aoaa iaprovaaant was notad in tha GI8S aodal (Soaarvilla at 
al., 1974) tha axact aituation is still pussling. Vhy should a change 
in raaolution at acalaa far reaoved froa tha bulk of K E produca tha 
partial iaprovaaant (Manaba, 1970)? 

A corollary to tha vaak K E haa baan its diffuaa distribution, 
figs. 4 and 5 show observation, tha GISS nodal distribution, and tha 
GLAS nodal distribution of K^, plottad as a sonal avsraga varsus praaaura 
and latitude. Evan nora striking is tha conversion of K E to K N which 
haa a characteristic dipole pattern in tha observations and lacks it in 
the previous nodals (Pig. 5). Corresponding results vara notad for tha 
National Canter for Atnospheric Research (1TCAR) nodal (Baker at al., 
1977, Pig. 9)* Tha changa in Pigs. 4 and 5 is substantial, sines tha 
negative portion of tha dipole is now batter located and of appropriate 
intensity. Sons conbination of tha radiative iaprovaaant and tha 
decreased saoothing haa produced a significant changa in those 
distributions, both in the integrated energetics and aaongst tha 


conversions. 


5* Spatial energetics 


Our cantral results concern the interplay of the eddy kinetic 
energy and its spectral subcomponents* Fig* 6 shows the observational 
and model values of K M and K E for the first two weeks of each time 
period. For two relatively different synoptic eituationa the model 
tracks the atmosphere's K E out to a period approaching 12 d. The 
agreement is not ao much one of quantitative precision aa trend 
following with a probability that is visibly non-random. Nc .e that for 
the Kg rise in January 1977 there appears to be a 1 d lag. A reasonable 
interpretation la the existence of predictive skill in an averaged 
quantity out to the predictability limits. 

While these two casss . eem stronger than a random correlation, the 
third canr uhows far less success. For January 1976 the model aeema 
unable to maintain or generate the Kg needed to follow the atmosphere' a 
rlae. Upon examining the energy flows more closely, we discover several 
significant phenomena. Consider firat tropospheric K*|. The cases 
consist of one low Km state (January 1976), one average Km state 
(January 1975), and one strikingly high Km state (January 1977). The 
latter case is unusual and was sustained at or near this high value 
throughout the 2 weeks. The model, in contrast, for the three cases has 
only one pattern t nonotonic growth. Some of the success of January 
1977 is clearly due to the atnoaphers's happening to match this built in 
bias of tho model. 
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An even stronger anomaly appear* In tha time displays of 
stratospheric K M shown in Pig. 7. There appears to be a ubiquitous 
tendency for sonal energy to collect in the model's topmost layer, with 
significant excess energy present in the model's one stratospheric level 
in periods of 4 to 6 d. Since there are no significant diasipativs 
mechanisms at that level, the excess energy appears to remain trapped 
there, functioning as a reaevoir of K M and K £ , but distorting the 
behavior of the individual wavenumbers. In another paper (Kalnay-Rivaa 
and Tenenbaum, I960) we examine this effect over the first 46 h end show 
the existence of major oscillations with a period of about 24 h. As 
discussed there, these oscillations do sesm to be attenuated by enhanced 
vertical resolution near the tropopause. 

The monotonic rise to unrealistically high levels also appears in 
stratospheric values of Kg (Pig. 6). The 1976 case gaint eddy kinetic 
energy in the model stratosphere much faster than the other two years, 
and at the same tine that the troposphere is rapidly losing energy. The 
failure in this csss is due to the model's inability to correctly 
maintain a strong polar jet stream over Norway, resulting perhaps from 
the smoothing at higher latitudes. 

Perhaps our most suggestive result concerns the longitudinal 

dependence of the stratospheric leaksga. In all three cases we see what 

appears to be vertical transport of ensrgy upward from the initial jets. 

Pig. 9 shows a time series of meridional cross sections of the sonal 
o 

wind at 90 E, the entry region of the subtropical jet centered south of 
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Japan* The obaarved data ahow noderate fluctuations while the nodal 
ahowa a converalon fron an "o" pattern to a "u" pattern with a 
characteriatic tine of 4 to 8 d. Thla aane behavior occurred above all 
Jets for all three caaea and aeena to be the antecedent of the ":i" 
ahaped average pattern characteriatic of noat nodela. (Fig. 4 above; 
Miyakoda at al., 1972, Fig. 2). 

An alternate way of exanlning thla phenonena ia ahown in Fig. 10 

which preaenta the tine aeriea of the aonal wind at the atratoapheric 

level. The obaervational aeriea ahowa noderate fluctuationa of regiona 

with wind a greater than 30 n/a centered above Norway and Japan. No 

valuea exceed 40 m/a. In atrong contrast, the nodel'a tine aeriea ahowa 
% 

an inkblot-like riae above Japan which growa to apeeda exceeding 60 n/a 
and which propagatea downstream to the Aleutle 3 (9 January) and the 
Yukon (13 January). The vertical leakage above ic atrong initial Jets 
occurred in all three yeara though only in the case illustrated (January 
1975) did the energy propagate downstream. For January 1976 and 1977 
the pool of energy renained above its origin (Norway and Japan, 
respectively). 

Thla behavior aeena very striking, and nay provide a clue to the 
cold polar atratoaphere phenomena seen in noat general circulation aod*l 
winter ainulatlona. Excess kinetic energy ia reaching the nid-latitude 
atratoaphere, and consiatent with the themal wind relation, erroneously 
modifying the polar stratosphere. 
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6. Spectral energetics 


In striking contrast to the successful K £ predictions for January 
1975 and January 1977 is the behavior of wavenumber 3. One can isolate 
the one or two dominant wavenumbers with little difficulty. Fig. 11 
shows the time history of Kj and K 3 for these two years. (These two 
wavenumbers are 100£ higher than the peak value of all other 
wavenumbers). In spite of the noisy signal ( one sees a much more 
discrepant behavior with gross errors in wavenumber 3 at tines when the 
summed Kg's remain correct. 

In terns of Saltzman's spectral energetics analysis one can imagine 
two possibilities: (l) erroneous transports in or out of the Northern 

Hemisphere troposphere, or (2) erroneous conversions from zonal 
potential or kinetic energy. As is inevitable, the sorting out of cause 
and effect will probably be difficult. 

As shown in Fig. 10 the erroneous stratospheric energy seems to be 
billowing up above the quasi-permanent jets. This transport shows up 
very strongly in both K M (Fig. 7) and Kg (Fig. 8). As a result 
wavenumber 3 has been deprived of the major portion of the energy needed 
to maintain the observed values of the subtropical jets in their 
characteristic standing wave pattern. Vavenunber 1 , which is dominated 
by a combination of the polar asymmetry and hemispheric contrasts, is 
less affected. 
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The spectral analysis approach is useful in pinpointing the 
location of other difficulties. As shown in Fig. 11 , the model's K 3 for 
January 1977 does very poorly starting about 6 January of the run. 
Examination of the various conversions away from the locations of the 
leakage to the stratosphere shows a clear example of how the spectral 
conversion can pin down the cause of this discrepancy in a causal sense. 
Fig. 12 shows the latitude-pressure distribution of K 3 , and the 
converaions C(K 3 ,K M ) and C(K 3 ,K n ) [Saltsman's M and L, respectively] for 
OOOO GMT on 1 January 1977 and 9 January 1977, Just before the start of 
the discrepancy. While the model has correctly predicted the 
qualitative and rough quantitative behavior of K 3 , the conversions are 
already qualitatively incorrect. C(K 3 ,K m ) is erroneously draining the 
model's K south of the sub-tropical Jet while C(K 3 ,K n ) is building K 3 
at what will be the location of the statoapheric warming 4 d later. The 
model displaces this conversion to coincide with the Jet itself. 

One other commonly suggested difficulty does not appear. The 
phases of the energy in individual wavenumbers appear to agree when the 
wavenumber magnitudes themselves agree. Fig. 13 shows the phase 
diagrams corresponding to Fig. 6 for January 1977. Both model and 
observation appear matched in phase until 7 January 

The spectral analysis techniques seem useful in on other way. 

Fifty percent of the amplitude of the stratospheric wanning of 11 
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January appears in observational wavenumber 3 analyses at approximately 
70°N. The model missed this warming, though in fairness, this was a 
feature not appearing until 10 d into the forecast. The presence of a 
very strong observational signal in wavenumber 3, and dominantly there, 
suggests a potential analytic tool for following up this phenomena. 

Two other general conclusions emerge from the spectral results. 

The first, not previously noted, concerns the behavior of the individual 
wavenumber Kj/a. One appears to have a phenomena best described as 
single wavenumber dominance. In Fig. 14 we show 30 d series of the 
first six observational K n . On a purely chance basis there appears to 
be far less than the expected amount of overlaps between the individual 
peaks. The relation of this result to some of Charney's work on 
oscillations between climate states is not clear. 

The second, noted by Tsay and Kao (1978) and others, is the very 
close tracking of the kinetic and potential eddy terms. As shown in 
Fi b , 15, the effect is clearly present for n • 1 and 3 and absent for n 
• 2. The other years show similar effects. It contrasts with the 
negative correlation shown by K M and K E in Fig. 6, a result which one 
tends to expect on the basis of energy conservation. 


7. Son* problem area* for apeotral energetics 


The Insights provided by the spectral energetics approach have been 
presented above as techniques for localising model failure at tinea when 
integrated eddy quantities are showing agreement. There remain, 
however, several problem areas which require further study to increase 
one's confidence in the analytic method, or the behavior of the 
atmosphere, or both. These areas include the issues of atmospheric 
cause and effect, the consistency of the conversions, and the 
calculation of vertical velocity dependent quantities. 

Tsay and Kao (1978) and Kao and Chi (1978) have made detailed 
studies for the period 1 December 1975 through 29 February 1976 in the 
wavenumber- time and wavenumber- frequency domains, respectively. 

Working with the observations, and predominantly at 500 mb between 30°N 
and 60°N, they attempt to trace the growth and decay of waves in the 
long (n ■ 1 to 5) and synoptic (n ■ 4 to 8) scales. They conclude that 
"the growth and decay of the kinetic energy of the long and synoptic 
scale waves are primarily controlled by the transport of kinetic energy 
to and from the wavea through non-linear wave interactions while the 
contribution to the kinetic energy conversion tends to balance the 
effects of the Reynolds and friotional stresses. " (Kao and Chi, 1978, 
abstract) . 

Their paper contains several well documented examples of the 
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correlation of changes in individual K n with the non-linear converaion 
tern. Ve question, however, whether such data demonstrates causation or 
ainply correlation. Their data tenda to show changes in conversions 
occurring at similar times as the change in the corresponding energy, 
rather than clearly preceding it as in Fig. 12 above. 

Beyond this issue, more general questions remain concerning the 
consistency of results obtained by applying Saltsman's spectral 
energetics analysis. The equations themselves are not in question but 
rather the quality of results from their application to either model or 
observational data. Saltaman'a equations provide a closed system with 
respect to conversion* into and out of individual K n . The algebriac sum 
of the conversions integrated over an appropriate space and time period 
should equal the aecular change in the corresponding K n . While such a 
result requires assuming the validity of dissipation calculations, D n , 
the lack of agreement gives us some sense of our overall confidence in 
the results. Because of the unavailabity of independent data for D n , 
Saltsman had to calculate it only as a residual. 

The results are not encouraging. Table 2 presents this type of 
calculation for a sampling of model and obsehrational runs. For most 
cases the sign of the change ia incorrect and in some coses the 
individual K n should have been driven solidly and unphysically negative. 
Three poaible explanations seem most likely: errors in C n , D n , and 

vertical transports. The conversion C n , and to a lesser extent some 
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others, depend on the vertical velocity, u . For the observations this 
is s notoriously difficult qusntity to evaluate. D n is s parsmeterised 
quantity whose physical basis is known in broad terns but whose detailed 
quantitative evaluation requires confidence in surface winds and 
boundary layer effects. Vertical transports are clearly a problem for 
the model and may be a problem for the observations. 

Ve will examine the vertical velocity problem in detail in a 
subsequent paper (Tenenbaum and Reddi, to be published). Our 
preliminary results show that alternate formulations for u (changing 
from mass convergence to iterative solution of the diagnostic omega 
equation can yield major changes in long- and synoptic scale values of 
C and extensive smoothing of the higher wavenumbers. 


8. Conclusions and future directions 


Ve have studied the model's medium range skill and climatology 
using various possible integrated quantities. Our major conclusions are 
as follows: 

(1) The model can show medium range skill in an integrated 
quantity, K , the Northern Hemisphere tropospheric eddy kinetic energy, 
out to the predicability limit of two weeks in two of three randomly 
chosen cases. 

(2) The model climatology is significantly Improved over previous 
versions in terms of the monthly time-mean energetics. This effect is 
most striking as sharper gradients in the sonal average of K E and the 
characteristic dipole behavior of the barotropic conversion terms. 

(?) The model retains a bias in accumulating energy in the topmost, 
stratospheric, level. The accumulation appears predominantly above the 
strong initial jets and downstream of the location of the 
quasi-permanent Jets associated with steep orography. The tropospheric 
K m also appears to grow to unreasonably large values. 

( 4 ) In the presence of this stratospheric leakage, wavenumber 3 
appears questionable in the role of an integrated energetics quantity: 
K e can show agreement while wavenumber 3 is erroneous. An expanded 
version of Ssltsman's spectral analysis (including transports) does 
provide a useful tool in showing cause and effect relations for where 
the model is failing to match significant energy flow in the atmosphere. 
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Then* approaches must be used with care aince Tsay and Kao 'a studies 
sees to show more correlation than causation. The spectral coefficients 
do show interesting results in terms of single wavenumber dominance (not 
previously noted) and kinetic energy tracking of potential energy (noted 
by Tsay and Kao). 

A number of future directions sre indicated. More synoptic cases 
should be studied, both randomly chosen, and specially chosen to explore 
the factors cited above (Km high or low, fast or slow propagation to the 
stratosphere, jet strong and well or poorly defined). A model with 
better resolution near the tropopause is needed which probably implies a 
shift to uneven spacing in pressure. If the stratospheric leakage is 
the result of erroneous waves induced by steep terrain, we may need to 
improve the conservation properties of the model in these regions (see 
Arakawa and Lamb, 1977). Alternatively, if the problem is erroneous jet 
stream or vertical velocity initialisation, we may need to examine more 
complex initialisation schemes to reduce the initial shock. In view of 
the presence of strong oscillations and strong leakage in all three 
casea, both effects are most likely present. 
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Table Captiona 


Table 1 


Table 2 


I 

I 


K , K , and K_ at 0000 GMT on indicated datea for obaervatione 
n n b 

and model predictions. Data are for the Northern Hemisphere 
troposphere (120 mb to 1000 mb). Units: 10^ J / m^. 

Changes in for observations and model predictions. Upper 
section shoes secular change over the 2 day period indicated. 
Lover section gives corresponding change in implied by the 
conversions averaged over the corresponding 2 day period. The 
first two conversions feed energy to Ky the laat two drain 
energy from Ky These results are quite typical. Data is for 
the Northern Hemisphere troposphere (120 mb to 100C mb). 
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Table 


conversion balance 



observation 

Model 

s 

(10 5 J / « 2 ) 

do 5 j / m 2 ) 

OOOO OUT 8 Jan. 1977 

1.5 

1.5 

0000 OUT 6 Jan. 1977 

0.6 

1.9 

chance 

♦0.0 

0.0 

average conversions 

(V / m 2 ) 

(V / » 2 ) 

N t0 S 

-0.? 

0.0 

*. «0K, 

0.0 

-0.? 


0.1 

0.1 

S to ^ 

0.2 

0.4 

net conversion 

-0.6 

-0.8 


do 5 j / « 2 > 

(10* 1 J / a 2 ; 

implied chance 

-1.0 

1 

-14 
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Fleur* Caption* 


FI*. 1 


Fi*. 2 


FI*, 9 


FI*. 4 


FI*. 


FI*. 6 


Constant praaaura chart at 200 mb for 0000 QNT 1 January 1977. 
Contour Interval la 120 m. 

Constant praaaura chart at 200 mb for 0000 GNT 7 January 1977. 
Contour Interval la 120 a. 

Constant preaaur* chart at J00 mb for 0000 GMT 11 January 1977. 
Contour Interval Is 120 m. 

Monthly mean of eddy kinetic ener*y, K . (a) Climatolo*lcal 

P* 

avara** for January from Pelxoto and Oort ( 1 0*^4 > . (b) GT3S 

model from Tenanbaum (1976). (c) NMC observations for January 

1079. (d) GLAS model prediction for January 1979. 

Monthly mean of the conversion K^, to K^. (a) Climatological 

average for January from Palxoto and Oort (1974). (b) GISS 

modal fro* Tanenbau* (1976). (c) NMC observations for January 

1979. (d) GLAS model prediction for January 1079 . 

Tima history of tropospheric mean kinetic energy. K , and 

n 

eddy kinetic energy, for the first two weeks of eaoh caae. 
(a) January 1979. (b) January 1076. (c) January 1077. 

90 


Units: 


S 2 

1CT J / m . 


Fig. 7. Time history of mesn kinetic energy, K , for the stratospheric 

n 

level (10 mb to 120 mb). Model data for all three cases and 
observational data for January 1975. The observational data 
for January 1976 and 1977 were similar to 1975. Units: 10^ J 

/ m^ bar. 

Fig. 8 Time history of eddy kinetic energy, K , for the stratospheric 

i E 

level (10 mb to 120 mb). Model data for all three cases and 
observational data for January 1975. The observational data 
for January 1976 and 1977 were similar to 1975. Units: 10^ J 

/ m^ bar. 


Fig. 9. Jet stream cross sections of the sonal wind at longitude 90°E . 
(a) Observations, (b) Model predictions. Data are for 0000 
GMT on 1 , 5, 9, and 13 January 1975, respectively. Regions 
above 30 m / s are shaded. 


Fig. 10. Surface display of the tonal wind at the stratospheric level 
(10 mb to 120 mb), (a) Observations, (b) Model predictions. 
Data are for 0000 GMT on 1, 5, 9, and 13 January 1975, 
respectively. Regions above 30 m / s are shaded. 

Fig. 11. Time history of eddy kinetic energy of wavenumbers 1 and 3 

0 
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Fig. 


Fig. 


Fig. 


Fig. 


(a) January 1975. (b) January 1977. Unita: 10^ J / m^. 

12. Pressure- latitude display of K^ t the conversion to K^, 
and the conversion K to K,. (s) The initial state at <XX)0 

n " 

GMT 1 January 1975. (b) The situation at 0000 GMT 5 January 

1975 Just before starts to diverge between model and 
observation, (c) The model prediction at 0000 GMT 5 Janaury 
1975. 

15. Polar diagram of magnitude and phase for observations 

and model. Data are the layer average from 220 mb to 550 mb 

(level 5 of the GLAS model) at 50 °N. This combination of 

pressure and latitude is near the peak of the distribution. 

Labels indicate days after initial state, (a) January 1975. 

5 2 

(b) January 1977. Units (radial): 10 J / m . 

14. Time history of the individual K n for n • 1 to 6 for Janaury 

1975. (a) Observations, (b) Model. Units: 10^ J / m^. 

15. Time history of the available potential energy, P n , and the 

eddy kinetic energy, K , for n ■ 1 to 5 for observational and 

n 

c ,2 

model data for Janaury 1975. Units: 10 J / m , 
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